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ABSTRACT 

We present predictions for the abundance of Ly-a emitters in hierarchical structure 
formation models. We use the GALFORM semi-analytical model of galaxy formation to 
explore the impact on the predicted counts of varying assumptions about the escape 
fraction of Ly-a photons, the redshift at which the universe reionized and the cosmo- 
logical density parameter. A model with a fixed escape fraction gives a remarkably 
good match to the observed counts over a wide redshift interval. The counts at bright 
fluxes are dominated by ongoing starbursts. We present predictions for the expected 
counts in a typical observation with the Multi-Unit Spectroscopic Explorer (MUSE) 
instrument proposed for the Very Large Telescope. 
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1 INTRODUCTION 

Dedicated narrow-band searches for Ly-a emitters have 
proven to be very efficient at detecting high-redshift galax- 
ies (e.g. Hu & McMahon et al. 2003). Objects found by this 
technique have to be confirmed spectroscopically, to rule out 
possible low redshift interlopers that may arise due to emis- 
sion lines other than Ly-a falling within the targetted wave- 
length interval. Nevertheless, a significant fraction of the 
detections appear to be bona-fide Ly-a emitters, and the 
number of objects accumulated to date by this technique in 
the redshift interval 2.4 < z < 6.6 is quite impressive. The 
Ly-a emission line is also found in a significant fraction of 
Lyman break galaxies (e.g. which are selected on the basis 
of their continuum emission. 

The ubiquity of the Ly-a line is at face value surpris- 
ing, given that it is resonantly scattered by atomic hydrogen, 
and so is easily absorbed by even a small amount of dust in 
a neutral gaseous medium. It is suspected that most Ly- 
a emitters have galactic winds (as is the case with Lyman 
break galaxies) which allow Ly-a photons to escape from 
the galaxy after only a limited number of resonant scatter- 
ings (Kunth et al. 1998; Pettini et al. 2001). The Ly-a line 
typically shows an asymmetric profile characteristic of such 
a process (e.g. Ahn 2004). The physics of this phenomenon 
is complicated, however, and remains poorly understood. 
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We present here the first predictions for the abundance 
of Ly-a emitters at different redshifts made using a model 
which follows the formation and evolution of galaxies in a 
hierarchical universe. In previous work, simple, ad-hoc pre- 
scriptions have been used to assign star formation rates to 
dark matter haloes (Haiman & Spaans 1999; Santos et al. 
2004). In this Letter, we use a semi-analytical model to make 
an ab initio calculation of the distribution of galaxy masses 
and star formation rates at different redshifts (e.g. Cole et al. 
1994; Kauffman et al. 1994; Baugh et al. 1998; Somerville & 
Primack 1999; Hatton et al. 2003). The model we use is able 
to reproduce the observed properties of galaxies both locally 
and at high redshift (ICofe et al. 20001 |Baugh et al. 2004| . 
The abundance of Ly-a emitters is sensitive to the adopted 
cosmological model and to astrophysical phenomena, such 
as the fraction of Ly-a photons escaping from galaxies and 
the distribution of galactic dust. Semi-analytical models are 
ideally suited to the exploration of such a parameter space. 



The semi-analytical model is described in Section |5] In 
Section^ we first present a compilation of the available ob- 
servational data on the abundance of Lyman-a emitters at 
different redshifts, and then compare these data with the 
predictions from our models. Finally, we present our conclu- 
sions in Section 0] 
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MODEL 
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178 


•^rcion. 


/esc. 


MUSE COUNTS 
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0.3 


0.7 


0.93 


10 


0.02 


70 


B 


0.3 


0.7 


0.93 


10 


0.1 


248 


C 


0.3 


0.7 


0.93 


10 


dust 


366 


D 


0.3 


0.7 


0.93 


6 


0.02 


163 


E 


Q.3 


0.7 


0.93 


20 


0.02 


77 


F 


0.2 


0.8 


1.15 


10 


0.02 


145 


G 


1 





0.52 


10 


0.02 


59 



Table 1. The parameters of the semi-analytical models for which 
the abundances of Ly-a emitters are predicted. The first column 
gives the model label. The next three columns give the basic cos- 
mological parameters: the density parameter, f2, the cosmological 
constant, A, and the amplitude of density perturbations, as spec- 
ified by trg (<T8 values are taken from Eke et al. 1996). In each 
case the baryon density is Qj, = 0.04 and the Hubble constant 
is Hq = 70kms — ^^Mpc -1 . Column five gives the redshift z r cion 
at which the universe is assumed to reionize. Column six gives 
the fraction / csc of Lyman-a photons that escape from the model 
galaxies. In the case of model C, the escape fraction is computed 
from a dust extinction model, as described in the text. The fi- 
nal column gives the number counts of Lyman-a emitters for a 
reference MUSE observation. This is the number of emitters in 1 
square arcminute in the redshift interval 2.8 < z < 6.7 brighter 
than 3.9 X 10~ 19 ergs~ 1 cm~ 2 . 



2 THE MODEL 

We use the semi-analytical model of galaxy formation, 
GALF0RM, to make predictions for the abundance of Ly-a 
emitters as a function of Ly-a flux and redshift. The GALF0RM 
model is described in full in Cole et al. (2000) and Benson et 
al. (2003); further details of the model used in this Letter are 
given in Baugh et al. (2004). The GALF0RM model computes 
the star formation histories for the whole galaxy population. 
The following steps are taken to compute the Ly-a emission 
from a model galaxy: (i) The number of Lyman continuum 
photons is computed from the star formation history in the 
model galaxy and the stellar initial mass function (IMF). In 
the Baugh et al. (2004) model, quiescent star formation in 
galactic disks produces stars with a Kennicutt (1998) IMF, 
whereas bursts of star formation triggered by galaxy mergers 
produce a flat ("top-heavy") IMF. (ii) The luminosity of the 
Ly-a line is computed assuming that all Lyman continuum 
photons are absorbed in HII regions and produce Ly-a pho- 
tons according to case B recombination (Osterbrock 1989). 
(iii) The observed Ly-a line emission depends on how many 
Ly-a photons escape from the galaxy. We have taken two 
approaches to estimating the escape fraction. In the first, 
we simply assume that a fixed fraction, / csc , of Ly-a pho- 
tons escape from the galaxy. Physically, this might arise if 
a fraction of the Ly-a photons escape through holes in the 
galactic gas and dust distribution. In the second approach, 
we calculate the absorption of Ly-a photons by a diffuse dust 
medium having the same spatial distribution as the stars, 
ignoring resonant scattering by neutral hydrogen. This may 
mimic what occurs if resonant scattering is suppressed by a 
galactic wind. We make a self-consistent calculation of the 
dust optical depth of the model galaxies, using the predicted 
galaxy scale length, gas mass and metallicity. Full details of 
this dust extinction model can be found in Cole et al. (2000). 
We plan a more detailed calculation of the escape fraction 
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Figure 1. The predicted number of Ly-a emitters in model A, 
compared to the observational data compiled in Table |2] The 
observational data are divided into different redshift ranges, as 
indicated by the lower key, and we have plotted only datapoints 
based on more than one galaxy. Model predictions are shown by 
lines, as indicated by the upper key. The thickness of the lines 
increases with redshift. The vertical and horizontal dashed lines 
indicate the sensitivity limits in flux and number density for a 
reference MUSE observation (see text). 

of Ly-a photons, including the effects of resonant scattering 
and gas outflows, in a future paper. 

We explore the impact on the abundances of Ly-a emit- 
ters of varying the redshift at which the Universe was reion- 
ized, Zreion, which is still poorly constrained. 

We examine the consequences of making three choices: 
(i) Zrcion = 10, (ii) Zrcion = 20 and (iii) « fe iot> = 6. The first 
two values are consistent with the optical depth to last scat- 
tering suggested by the WMAP measurement of the corre- 
lation between microwave background temperature and po- 
larization | |Kogut et al. 2003t . The latter value is suggested 
by the detection of a Gunn- Peterson trough in the spectrum 
of a z — 6.28 quasar by Becker et al. (2001). In the model, 
gas is prevented from cooling in haloes with circular veloci- 
ties below GOkms^ 1 for redshifts z < z r cion- The models that 
we consider are listed in Table Q Models A-E reproduce the 
observed luminosity function for the local galaxy population 
(e.g. Baugh et al. 2004). Finally, we also consider the effect 
on the counts of Ly-a emitters of varying the cosmological 
density parameter, £1 (models F and G). Note that we have 
not attempted to vary any other parameters of the GALF0RM 
model in these two cases to force the model to reproduce 
the local galaxy luminosity function. 



3 THE OBSERVATIONAL DATA AND MODEL 
PREDICTIONS 

We list in Table [2] a compilation of published observational 
data on number counts of Ly-a emitters at different observed 
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Table 2. Data Compilation. The data are divided into unit redshift intervals: the following symbols are used to denote data from each 
redshift interval in the figures (■ :[z < 3], T :[3 < z < 4], A :[4 < z < 5], : [ z = 5.1], «:[5 < z < 6], X:[6 < z < 7]). Col.l: redshift; 
Col. 2: redshift interval; Col. 3: Ly-a flux (in 10~ 17 ergscm~ 2 s~ 1 ); Col. 4: cumulative counts per unit solid angle per unit redshift (in 
arcmin" 2 ); Col. 5: Poisson error on counts (in arcmin -2 ); Col. 6: number of Ly-a emitters; Col. 7: area of survey (in arcmin 2 ); Col. 8: 
factor applied to correct for contamination by low-z interlopers; Col. 9: method (NBF=narrow band filter, LS=long-slit spectroscopy); 
Col. 10: method used to reject or correct for low-z interlopers (EW=equivalent width, spec on N = follow-up spectroscopy of N objects); 
Col.ll: reference (DOl: Dawson et al. 2001; F03: Fujita et al. 2003; H98: Hu et al. 1998; H02: Hu et al. 2002; K03: Kodaira et al. 2003; 
KOO: Kudritzki et al 2000; O03: Ouchi et al. 2003; ROO: Rhoads et al. 2000; R03: Rhoads et al. 2003; Sa04: Santos et al. 2004; S03: 
Shimasaku et al. 2003; S04: Shimasaku et al. 2004; SteOO: Steidel et al. 2000; StiOl: Stiavelli et al. 2001) 

* i„A r„„ t — + — c — a ** correc ted for gravitational lensing 
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Figure 2. The predicted number counts of Ly-a emitters in 
model A, in the redshift range z = 3 — 4. The contributions from 
quiescent star forming galaxies and ongoing starbursts are plotted 
separately, as indicated in the key. The vertical and horizontal 
dashed lines are the same as in Fig. Q 



line fluxes and redshifts from "blank field" surveys that we 
will compare against our model predictions. Note that some 
of the quantities listed in Tableware derived from the infor- 
mation given in the original sources, so we present this com- 
pilation in order to facilitate future comparisons between 
data and model predictions by other authors. We do not 
include data from surveys targetted around high-z objects 
such as quasars (e.g. Fynbo et al. 2001), since in this case the 
number density of Ly-a emitters may be biased by an un- 
known factor. A few of the surveys used spectroscopy to di- 
rectly search for Ly-a emitters. However, most used narrow- 
band imaging to identify objects having a strong emission 
line at the wavelength corresponding to Ly-a at a partic- 
ular redshift. Samples of objects obtained in this way are 
generally contaminated by lower redshift galaxies for which 
some other emission line (e.g. [OJ7]3737, [OIII]5007, Ha) 
happens to fall at the same wavelength. This contamina- 
tion is typically estimated and removed using the equivalent 
width of the emission line and/or broad-band colours, or 
from follow-up spectroscopy on a subsample of the objects. 
The methods used in the different surveys are indicated in 
the table. 

The observational data are plotted in Fig.0 with differ- 
ent symbols indicating data from a given unit redshift range. 
The predictions of our fiducial model, A from Table are 
shown by the solid lines in Fig. The escape fraction / eac 
was set to give a reasonable match to the observed number 
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Figure 3. The predicted counts of Ly-a emitting galaxies, plotted in different redshift intervals, as indicated by the key. (a) Shows the 
predictions for model A; a subset of these are reproduced in each panel for reference. The data from Fig. 1 are also plotted here using the 
same symbols as before. The remaining panels show the impact on the predictions of changing different aspects of the model: (b) Varying 
the escape fraction or using a dust model, (c) Varying the redshift of reionization. (d) Varying the cosmological density parameter. 



counts at z ~ 3. This simple model, in which the escape frac- 
tion is independent of redshift and galaxy properties, does 
a surprisingly good job of matching the observed counts at 
different redshifts. 

The separate contributions to the counts of Ly-a emit- 
ters from galaxies that are forming stars quiescently and 
from ongoing starbursts (which in our model are triggered 
by galaxy mergers) are shown in Fig. [5]for model A over the 
redshift interval 2 = 3 — 4. Quiescently star forming galax- 
ies dominate at fainter fluxes, whereas starbursts account for 
the brighter Ly-a sources. This is largely due to the flat IMF 
assumed in starbursts, which typically yields ten times the 
number of Lyman continuum photons for a given amount 
of star formation, compared with the Kennicutt (1998) IMF 
that we adopt for quiescent star formation. 

In Fig. |3J we present the predicted counts of Ly-a emit- 
ters as a function of redshift for the different GALF0RM models 
listed in Tabled The results for model A, our fiducial model, 



are reproduced for reference in each panel. In Fig. |3Jb) we 
show the effect of increasing the escape fraction / osc by a 
factor of 5. Fig. O^b) also shows that assuming a fixed es- 
cape fraction / csc = 0.1 produces similar number counts to 
a model in which the Ly-a photons are absorbed by dif- 
fuse dust without resonant scattering, as might occur in a 
galactic wind. Fig. EJc) shows the impact on the predicted 
counts of varying the redshift at which the universe is reion- 
ized. For model D, with z rc ion = 6, the number of faint Ly-a 
sources changes substantially for 5 < z < 6 and 4 < z < 5, 
compared with the cases where z re ion ^ 10. This is because 
in model D gas is still able to cool in low mass haloes (i.e. 
with circular velocities below 60kms _1 ) up to z = 6, and is 
available to form stars and thus generate Lyman continuum 
photons. Finally, Fig. OJd) illustrates the effect of changing 
the cosmological density parameter whilst retaining a flat 
universe. The changes in the model predictions in this case 
primarily reflect the change in the normalisation of density 
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fluctuations, as specified by as in Table^ which is adjusted 
to reproduce the local abundance of rich clusters when the 
density parameter is varied. 



4 DISCUSSION AND CONCLUSIONS 

The conclusion of this first study to look at the predictions 
of hierarchical models for the number of Ly-a emitters is 
that simple models do remarkably well at reproducing the 
observed counts. Assuming that typically just 2 % of the 
Ly-a photons escape from high-z galaxies, a value chosen to 
match the observed counts at z ~ 3, is sufficient to give a 
reasonable match to the observed counts at faint fluxes over 
the redshift interval 2 < z < 6. 

This study demonstrates the capability of semi- 
analytical modelling to make predictions that can serve as 
an input into the design of new instruments. The Multi- 
Unit Spectroscopic Explorer (MUSE) (B acon et al. 20041 has 
been proposed to ESO as a second-generation instrument 
for the Very Large Telescope. MUSE will be able to identify 
Ly-a emitters over a redshift interval 2.8 < z < 6.7 over 
a field of view of f arcmin 2 . An exposure of 80 hours will 
reach a 5 a sensitivity of 3.9 x 10 -19 erg/s/cm 2 . We predict 
that MUSE will be able to detect a large number of such 
objects at this flux limit: around 70-400 per arcmin 2 (see 
the final column of Table 0. Observations with MUSE will be 
able to exclude some of the models we have considered, and 
therefore remove some of the uncertainties in our modelling 
of Ly-a emission. More importantly, such observations will 
provide a critical test of our ideas about star formation in 
objects at high redshifts. 
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